Reply to Comment by A. Moroz 


Sajeev John and Ovidiu Toader 
Department of Physics, University of Toronto, 
60 St. George Street, Toronto, Ontario, Canada M5S-1A7 
(January 19, 2000) 
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In his comment, Moroz questions the validity of the 
near band edge (effective mass) approximation to the to- 
tal photon density of states (DOS) as a useful represen- 
tation of the local density of states (LDOS) experienced 
by a single radiating atom or molecule located at a par- 
ticular position r within a photonic crystal (PC). In this 
approximation, the band edge DOS takes the form: 


const cj — lu c 


where 77 = —0.5 for a 1-d PC and 77 = 0.5 for a 3-d PC. 
We reassert that this behaviour indeed applies to the 
LDOS as well as the DOS. However, the frequency range 
over which this behaviour is realized depends sensitively 
on r. In particular, if r is chosen near a node of the 


E(r) 


then uj must be 


electromagnetic field intensity 

chosen very close to to c before the asymptotic behaviour 
is realized. The seemingly arbitrary exponents obtained 
by Moroz are simply an artifact of fitting the asymptotic 
form to numerical data for a frequency lo which is not 
sufficiently close to lu c at certain positions r. 

We consider precisely the example quoted by Moroz in 
his comment and assume that the LDOS has the asymp- 
totic form: 

p(w,r) = fC(r)\w c -u\ v 

Near the lower band edge of the first photonic band gap 
(lu ^ lo c ) we define u = 1 — — > 0. In order to numeri- 
cally estimate the exponent 77, we write: 


V = logio P = *?(log 10 u + log! 


k>gio K(r) 


Using equations (4) and (7) of Moroz's paper M we plot 
(in Fig. la) y as a function of z = log 10 u for 8 different 
positions r in the 1-d unit cell of the example quoted in 
the above comment. The asymptotic behaviour of dy/dz 
for large negative values of z (cu — > uj c ) yields the expo- 
nent 77 (see Fig. lb). In this model the lower band edge 
mode intensity vanishes at x = \f\ = 0.5 (center of air re- 
gion) and has a maximum at x = 0.0 (center of dielectric 
slab). For all cases the asymptotic behaviour (u> — » u) c ) 
yields the common exponent 77 = —0.5. However arbi- 
trary values of dy/dz, and hence 77, may be erroneously 
inferred by choosing too large a value of \ui — oj c \. This is 
particularly evident near the node of the field intensity. 

We conclude that although the LDOS is sensitive to 
the actual position the exponent 77 is indeed universal 
except on a set of measure zero, namely the field inten- 
sity nodes. The seemingly arbitrary exponents quoted 



by Moroz are somewhat misleading. On the other hand, 
inhomogeneous line broadening is a very important and 
relevant ingredient which must be incorporated into the- 
oretical models which aim to interpret experiments in- 
volving a distribution of atoms in a PC. 
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FIG. 1. 


z=log ]() u 

Top picture shows log 10 p (u>, x) as a function of 
z = log 10 u for 8 positions in the unit cell. Here x — \f\. 
Bottom picture shows the slope of the curves. In all cases the 
asymptotic behaviour (ui — ■> ui c ) yields r\ — —0.5. 
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